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I.  EXPERIMENTAL  WORK 

In  the  Final  Report  for  1974  we  described  at  length  sogie  detailed 
results  on  relaxation  times  In  monolayer  and  sub-monolayer  He  films 
adsorbed  on  Vycor  porous  glass,  together  with  the  situation  concerning 
experiments  on  the  density  of  adsorbed  helium  as  deduced  from  capacitance 
measurements  and  progress  on  the  modified  Pomeranchuk  cooling  project. 

All  three  sets  of  experiments  have  been  continued  with  considerable  success, 
and  we  present  here  a summary  of  results  and  of  analysis  as  so  far  completed. 

1.1  Dielectric  constant  of  adsorbed  helium 


In  previous  reports,  we  have  described  this  work  as  a measurement 
of  the  density  of  adsorbed  helium,  but  the  results  presented  below  lead  to  the 
conclusion  that  the  correlation  between  density  and  dielectric  constant  is  not 
understood,  and  we  therefore  now  prefer  to  describe  them  only  as  dielectric 
constant  experiments. 

As  already  explained,  the  experimental  arrangement  is  quite  simple, 
and  consists  essentially  of  a parallel-plate  capacitor  containing  Vycor  porous 
glass  into  which  helium  can  be  condensed  at  pressures  between  0 and  about 
35  atm.,  and  temperatures  between  IK  and  4K.  The  difficulties  which  we 
described  in  Section  1,1  of  the  last  report  have  been  eliminated  by  the  means 
of  shielding  proposed  there,  and  we  have  confidence  in  the  validity  of  our 
results.  Figure  1 is  a diagram  of  the  final  design  of  the  low  temperature 
capacitance  cell.  Measurements  of  capacitance  were  made  with  a MOSFET 
oscillator  at  a frequency  20MHz.  The  effective  dielectric  constant  is  the 
ratio  of  the  capacitances  with  and  without  dielectric  material  between  the  plates, 
and  because  of  the  porous  nature  of  the  material  this  has  then  to  be  Interpreted 
to  give  the  dielectric  constant  of  the  body  of  the  porous  material,  and  of  the 
helium  in  the  pores . 


Figure  2 is  a plot  of  the  dielectric  constant  C ' of  porous  Vycor  and 
of  c v of  the  Vycor  glass  material.  The  first  of  these  is  derived  directly  from 
the  measurements,  the  second  is  given  by 
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where  uis  the  fractional  volume  of  the  voids,  i.e.  the  porosity, 
stated,  the  temperature  variation  is  unusually  large,  for  a reason  which  we 
cannot  explain.  We  are  at  present  examining  various  possibilities. 


(1) 

As  already 


4 

A He  adsorption  isotherm  whl^h  was  carried  out  ■ T=1.995K  is  shown 
in  Figure  (13).  Equal  quantities  of  the  He  were  condensed  Into  the  cell,  each 
quantity  being  represented  by  a division  on  the  abscissae.  They  condensed 
from  a fixed  volume  of  pipework  in  a period  determined  by  the  time  taken  for 
the  pressure,  measured  with  an  oil  manometer,  to  drop  between  fixed  limits. 
The  vapour  pressure  of  the  sample  was  measured  with  a Wallace  and  Tiernan 
Gauge  which  was  calibrated  against  a mercury  manometer  and  found  to  be 
accurate  U>  within  } torr.  The  ’dielectric  constant’  values  were  calculated 
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using  a value  of  u equal  to  that  used  in  analysing  the  fille|  pores  data  (see 
below);  thus  (apart  from  the  points  which  correspond  to  He  with  filled  pores) 
the  ordinate  does  not  represent  a true  dielectric  constant.  The  change  in 
slope  which  occurs  In  the  "dielectric  constant"  curve  is  characteristic  and 
is  probably  associated  with  the  completion  of  the  first  layer. 

4 

Values  of  dielectric  constant  for  He  "filled  pores"  are  shown  in 
Figure  (4).  The  value  of  u wa|  taken  to  be  0.3  from  our  previous  work. 

Suitably  normalized,  our  bulk  He  data  are  also  shown  (some  of  the  points 
taken  around  T.  are  omitted  for  clarity)  and  are  compared  with  the  values 
of  the  dielectric  constant  obtained  from  the  known  density  measurements. 

The  well  established  value  for  the  polarizability,  0.123  (c.g.s.  units)  was 
used  to  convert  the  density.  Agreement  of  the  bulk  values  Is  extremely  good, 
and  gives  us  confidence  in  our  basic  results.  The  points  for  helium  in  Vycor 
show  an  increasing  discrepancy  as  the  temperature  increases,  which  must  be 
attributed  to  the  characteristic  properties  of  the  restricted  geometry.  The 
density  maximum  is  much  closer  to  the  normal  lambda  temperature  than  are 
the  temperature  of  the  specific  heat  maximum,  or  that  of  onset  of  superfluidity. 

4 

The  dielectric  constant  was  also  measured  for  He  at  near  monolayer 
coverage.  The  quantity  of  4 He  to  condense  was  estimated  from  the  adsorption 
isotherm  (Figure  (3),  and  corresponded  to  two  divisions  along  the  abscissae. 
This  may  be  less  than  the  quantity  needed  for  a complete  monolayer,  especially 
if  only  the  'dielectric  constant'  curve  is  taken  as  a guide  (see  above).  The 
analysis  required  an  estimate  of  u for  a 4He  monolayer.  By  assuming  It  to 
consist  of  a single  layer  of  atoms  3.6  A in  diameter,  and  using  the  value  for 
the  specific  surface  area  obtained  previously  in  this  laboratory,  a value  of 
0.06  was  deduced.  Values  of  the  dielectric  constant  obtained  using  this  figure, 
but  uncorrected  for  the  effect  of  desorption,  are  shown  in  Figure  (15).  This 
correction  was  important  because  the  dead  spaces,  especially  those  at  low 
temperatures,  were  relatively  large.  It  was  also  difficult  to  estimate  because 
of  the  complex  cell  geometry,  and  the  effect  of  the  changing  temperature  profile 
along  the  sample  lines  as  the  level  of  liquid  He  in  the  *He  dewar  dropped. 
Because  of  these  difficulties  the  correction  was  only  accurate  enough  to  be 
indicative.  It  was  only  worked  out  at  4K,  and  is  represented  by  the  arrow  in 
Figure  (6);  it  is  evident  that  the  temperature  dependence  for  a *He  monolayer 
is  at  least  small. 

4 

In  absolute  terms,  both  the  He  'full  pores'  and  monolayer  dielectric 
constants  are  interesting  because  they  are  both  higher  than  can  be  accounted 
for  by  the  increased  density  of  He  in  Vycor.  This  might  be  due  to  the  effects 
of  the  particular  geometry  of  the  Vycor  cavities,  but  the  measurements  under 
pressure  described  below  suggest  that  another  explanation  may  be  correct. 

Measure  me  qts  of  the  dielectric  constant  as  a function  of  pressure  were 
carried  out  at  the  temperatures  1.41K  and  2.08K.  The  pressure,  which  waB 
applied  via  ope  of  the  two  sample  lines,  was  measured  with  a Budenberg  dial 
gauge  connected  to  the  other  line.  Pressures  could  be  read  to  within  i p.s.l. 

In  this  way  it  was  possible  to  correlate  any  sudden  changes  in  frequency 
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(indicative  of  a He  liquid-solid  transition  in  Vycor)  with  blockage  of  the 
sample  lines  due  to  the  formation  of  solid  bulk  ^He.  The  results  are  shown 
in  Figure  6,  calculated  according  to  two  models,  in  one  of  which  ("single 
sphere")  the  filled  Vycor  pores  are  assumed  to  behave  as  a collection  of 
non-interacting  spheres,  while  the  other  assumes  that  the  internal  field 
connection  is  zero.  In  either  case,  the  variation  in  c is  much  less  than 
would  be  expected  from  the  statistical  layer  model  of  the  film,  unless  60% 
of  it  were  incompressible.  We  thus  need  to  seek  another  explanation. 


Another  possible  explanation  is  that  an  induced  dipole  exists  in  the 
atoms  of  the  first  monolayer,  resulting  from  the  very  strong  adhesive  forces, 
which  lead  to  the  enhanced  dielectric  constant  of  the  first  layer  and  the  filled 
pores.  A difficulty  with  this  suggestion  is  that  it  is  not  clear  in  detail  how 
the  induced  dipoles  could  orient  in  the  applied  eleotric  field,  since  they  may 
be  expected  to  remain  always  perpendicular  to  the  surface.  It  Is  possible  that 
this  could  occur  if  the  atoms  move  slightly  in  a sharply  varying  surface  potential. 
If  so,  then  the  anomalously  small  variation  of  dielectric  constant  with  pressure 
shown  in  Figure  6 must  be  due  to  inhibition  of  the  dipole  orientation  as  the 
pressure  increases.  This  somewhat  speculative  explanation  is  currently  being 
further  explored. 


1.2  Pomeranchuk  refrigeration  with  adsorbed  He 


In  the  past  year,  we  have  firstly  carried  out  with  complete  success 
conventional  Pomeranchuk  cooling  of  liquid  3He  into  the  superfluid  phases, 
using  a reconstructed  apparatus.  Secondly  we  have  installed  the  off-melting 
curve  cell  and  carried  out  experiments  which  appear  to  show  (t)  that  the  basic 
requirements  for  the  process  are  satisfied,  allowing  us  to  look  for  the  polar 
phase  of  superfluid  He  without  solid  being  present,  (li)  that  the  range  of 
Pomeranchuk  cooling  is  significantly  extended,  and  (ill)  that  the  method  is 
particularly  advantageous  for  adapting  to  indirect  cooling. 


For  the  first  of  these  experiments  we  made  several  small  but  important 
alternations  to  the  apparatus,  including  improved  thermal  anchoring  of  the  ^He 
compression  line  and  enlarging  the  channel  connecting  the  mai^i  compresslonal 
cell  to  the  strain  gauge,  which  had  become  blocked  with  solid  He.  These 
improved  the  performance  considerably,  but  we  were  still  unable  to  cool  below 
about  5mK.  After  several  further  experiments  which  included  examination  of 
the  operation  of  the  bellows  at  room  temperature,  it  became  clear  that  the 
particular  bellows  configuration  used  was  unstable.  In  order  to  reduce  as  far 
as  possible  the  dead  space  occupied  by  bulk  liquid  3 He,  they  had  been  designed 
so  that  compression  of  the  3He  was  achieved  by  compressing  the  bellows  themselves. 
In  this  case,  it  is  possible  for  them  to  move  sideways  and  scrape  against  the 
side  of  the  container  with  considerable  frictional  heating.  A more  stable  (and 
conventional)  arrangement  is  to  compress  the  3He  by  extending  the  bellows 
(l.e.  the  He  is  on  the  outside  rather  than  the  inside).  When  the  stable 
arrangement  was  used  the  main  frictional  heating  problem  disappeared.  In 
several  subsequent  experiments  we  were  able  to  reduce  the  temperature  at 
melting  pressures  from  about  32mK  to  below  1.2mK  (i.e.  through  both  the  A 
and  the  B transitions  of  superfluid  3He)  and  keep  it  there  for  periods  of  the 
order  of  an  hour.  The  conventional  Pomeranchuk  ooollng  has  therefore  been 
completely  successful. 


^ The  dielectric  constant  of  THE  He*  MONOL.AVER 
(UNCORRECTED  FOR  THE  EFFECT  OF  DESORPTION) 
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The  next  step  was  to  replace  the  conventional  cell  with  the  specially 
designed  version  containing  Vycor  porous  glass  discs  and  incorporating  a 
hydraulically  operated  needle  valve  on  the^He  cell.  With  this  apparatus  we 
were  able  to  carry  out  a series  of  experimen  ts  to  measure  the  pVT  relation 
of  the  combined  bulk  liquid/adsorbed  3He  system.  The  pressure  was  monitored 
in  the  usual  way  with  a strain  gauge  pressure  transducer,  and  volume  changes 
by  the  movement  of  the  bellows,  again  measured  with  a capacitive  technique. 
Solidification  in  the  Vycor  pores  should  be  detected  by  an  apparent  increase  in 
the  compressibility  at  pressures  below  the  melting  pressure.  Unfortunately 
during  these  runs  the  dilution  refrigerator  was  not  operating  properly,  and 
we  were  able  to  take  measurements  only  at  about  95mK.  These  did,  however, 
seem  to  confirm  that  there  is  a melting  line  for  the  adsorbed  3He  which  is  about 
10  p.s.i.  below  the  bulk  melting  curve  at  the  same  temperature. 

Subsequently  we  carried  out  a Pomeranchuk  experiment  which  gave  a 
cooling  of  15mK  without  formation  of  bulk  solid,  with  only  10T  of  the  total  3He 
in  the  Vycor.  For  various  reasons  we  were  again  not  able  to  start  compression 
at  a lower  temperature  than  105mK,  but  we  reached  a low  temperature  of  20mK, 
the  last  part  of  the  process  on  the  bulk  melting  curve.  This  represents  an 
important  enhancement  over  conventional  Poleranchuk  cooling  which  we  hope 
will  continue  into  the  lower  temperature  region  below  20mK.  The  reasons  for 
the  high  starting  temperature  are  not  entirely  clear  at  present. 

The  success  of  these  experiments  has  emphasized  an  important  point, 
namely  that  this  modified  cooling  process  provides  an  outstanding  advantage 
for  indirect  cooling,  and  we  shall  Incorporate  this  into  our  programme.  The 
advantage  arises  from  the  fact  that  the  inhibition  of  bulk  solid  formation  around 
the  heat  transfer  surfaces  prevents  the  consequent  reduction  in  interface 
conductance,  which  is  a serious  problem  in  the  conventional  method. 

The  results  of  the  pVT  measurements  outlined  above  were  presented 
at  the  Manchester  Solid  State  Physics  Conference  in  January  197o  (unpublished). 


1.3  Nuclear  magnetic  resonance  measurements 

We  have  continued  our  extensive  measurements  of  spin-spin  and  spin- 
lattice  relaxation  times  in  monolayer  and  submonolayer  3He  films  on  Vycor, 
and  further  analysed  our  previous  measurements  of  the  susceptibility  and  spin 
diffusion  coefficient.  These  are  reported  in  detail  in  Sections  III  and  II 
respectively. 
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n.  Some  Considerations  of  the  Magnetic  Properties  of  Multilayer  ^He  Films 

Adsorbed  on  Vycor  Glaas 
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Introduction 

The  following  notes  summarise  soveral  phenomenological  models  which 
have  been  developed  during  the  last  year  in  a n attempt  to  explain  some  of  the 
properties  of  liquid  ^He  contained  in  porous  Vycor  glass.  Multilayer  He  films 
are  fascinating  systems  in  which  layers  of  very  different  quantum  degeneracy 
exi9t  In  close  proximity.  For  example,  it  is  not  known,  when  inter-layer  motion 
occurs,  whether  particles  carry  their  moan  thermal  energy  and  'susceptibility' 
with  them  (quantum  assimilation)  or  even  whether  the  Landau  approach  can  be 
used  in  situations  whore  translational  invariance  is  violated^  K Effects  which 
may  be  relevant  to  superfluid  liquid  '*He  are  also  apparent. 

It  will  be  clear  from  these  notes  that  more  experimental  data  is  needed 
to  discriminate  between  the  varloa  s models  discufasedl  and  several  experiments 
which  may  prove  definitive  are  suggested. 

3 

I.  The  Magnetic  Properties  of  Multilayer  He  Films 

,2) 

Adsorption  isotherm  experiments  have  established  that  the  first  adsorbed 

3 

layer  has  a very  high  density  appropriate  to  solid  He  at  400  atm.  Thermal 

measurements^  suggest  that  the  dominant  excitations  are  phonons  and  N.M.R. 

measurements^  show  that  the  correlation  time  for  atomic  motion  within  the 

-7  -6 

first  layer  is  between  10  and  10  seconds,  both  pieces  of  evidence  endorsing  the 

conclusion  that  the  first  layer  is  solid  like.  In  view  of  the  tight-binding  to  the  wall 

(~50  K),  atoms  are  precluded  from  jumping  out  of  this  layer,  and  the  spatial  extent 

of  thoir  wavefunctious  is  probably  so  restricted  as  to  inhibit  significant  exchange 

coupling  with  second  layer  atoms.  The  properties  of  this  first  layer  are  believed, 

therefore,  to  be  substantially  unaffected  by  adsorption  of  a second  layer. 

(5) 

Magnetic  experiments  on  the  second  layer'  ' have  boon  carried  out  by  pre- 

4 

plating  with  a He  layer  and  it  is  found  the  atoms  are  much  more  mobile,  with  a 
-9 

correlation  time  ~10  seconds.  The  density  of  the  second  layer  is  equivalent  to 

(2) 

llquld/solld  on  the  melting  curve  (32  atm.)'  . The  susceptibility  of  a two-layer 

3Ile  film  is  best  explained^  by  assuming  the  first  layer  to  bo  magnetically 

independent,  and  treating  the  second  layer  as  an  antiferromagnetic  Curie-Weiss 

solid  with  Tvt  = 36  mdeg. , which  agrees  with  correlation  time  above  if  the  dominant 
N 


motion  is  exchange.  Adsorption  isotherm  measurements  have  shown  that  following 

3 

the  addition  of  more  layers  (which  are  essentially  like  bulk  liquid  Ho  under  the 

(2) 

S.V.P. ) the  second  layer  itself  is  depleted  somewhat'  ' and  may  therefore  become 

liquid-like  in  nature.  In  fact  down  to  0.4  K the  total  susceptibility  yean  be  adequately 

(6  7) 

described  by  the  expression  ' 

X Cl*l  + C2X2  + Vt  (1) 

where  c is  the  relevant  mass  in  each  region  of  the  film,  vis  the  (Curie)  susceptibility 

1 3 * 

of  the  monolayer,  v the  susceptibility  of  liquid  He  (extrapolated)  at  32  atm.  and 

r gv 

the  bulk  liquid  susceptibility  under  the  S.V.P. ' . Treating  tte  second  layer 
on  a molecular  field  modul  with  = 36  mdog.  of  coarse  works  equally  well. 

In  either  case,  to  get  perfect  agreement  (Fig.  1)  one  has  to  allow  for  a slight  density 
enhancement  of  the  third  layer;  the  magnetic  profile  is  still  changing  slightly  until  the 
fourth  layer  is  reached^’ 

Bolow  0.4  K,  however,  the  prediction  of  the  liquid-model  = bulk  liquid 
susceptibility  at  32  atm)  is  significantly  lower  than  the  experimental  value  (Fig.  1). 

(7) 

Three  explanations  have  boon  put  forward.  The  first'  attributes  the  discrepancy  to  a 

partly  solidified  second  layer  with  the  solidified  part  having  a Curie  susceptibility. 

Second,  the  second  layer  might  be  treated  ns  a homogeneous  Curle-Welss  solid, 

(6) 

although  at  the  lowest  tomperaturos  T ~T^  so  exact  calculations  are  difficult  \ 

Third,  the  second  layer  is  treated  as  a liquid  but  Is  believed  to  enhance  the  susceptibillty 
of  the  central  liquid  region  via  a substrate-mediated  Interaction^  In  section  2 
we  discuss  these  models  from  a different  standpoint,  with  particular  reference  to  the  last 
because  of  its  relevance  to  ultra-low  temperature  experiments  on  superfluid  1 He^1^. 

3 

The  properties  of  the  spin  relaxation  times  In  the  saturated  (fu!l-pores),  He 

film  are  Interesting  from  a number  of  points  of  view,  but  particularly  so  In  relation 

3 (11  12) 

to  the  behaviour  of  those  same  parameters  in  a single  Ho  monolayer'  * . Both 

3 rt  3\ 

T and  T,  are  very  much  shorter  than  the  bulk  liquid  He  values^  ' , and  arc,  In  fact, 

1 2 

comparable  with  the  values  seen  In  the  monolayer  (Figs.  2-5).  It  follows  that 
the  entire  film  is  relaxing  (or  dephasing)  as  a result  of  processes  taking  place  In  the 
adsorbed  layerl?\  The  Interesting  question  of  course,  is  how  the  spin  information 
is  communicated  to  the  rest  of  the  film  In  such  a way  that  the  magnetic  profile  Is 
preserved. 

A number  of  features  should  be  mentioned.  T^  for  the  monolayer  Is  very 
unusual  In  that  T^  Is  vaguoly  linear  In  frequency  (Fig  2),  also  passing  through  the  origin 


- 8 - 


and  therefore  making  It  difficult  to  estimate  the  correlation  time  (tc)  and  second  moment 
(Mg).  The  linear  behaviour  is  believed  to  bo  related  either  to  a range  of  T^s 
resulting  from  surface  heterogeneity*  *,  or  to  the  anisotropy  of  the  dipole  inter- 
action which  becomes  important  because  atoms  are  constrained  to  move  on  planes  at 
arbitrary  angles  to  the  main-field  direction*12, 15).  In  the  saturated  film  ^ goes 

o 

as  (frequency)  which  is  characteristic  of  most  liquid-like  (ut^  « 1)  systems 
(Fig.  3).  Since,  however,  the  monolayer  Is  still  the  origin  of  relaxation,  either 
the  heterogeneity  h;is  been  smoothed  out  or  the  dimensional  constraint  has  been 
removed  and  the  system  is  behaving  as  a truly  three  dimensional  entity.  ^ recoveries 
are  slightly  non-exponential*  ’ * Is  independent  of  temperature  below  1 K*  *. 

In  the  monolayer  T„  rises  ("10/3"  fashion;  as  frequency  is  increased  (Fig.  4) 

mt 

However  the  change  takes  place  in  a regime  where  T^»  Tg  so  that  it  is  difficult 

to  attribute  the  effect  to  life-time  broadening*12*.  In  the  saturated  film  (Fig.  5) 

a single  well-defined  T exists  but  there  is  no  such  increase  with  frequency,  and 

in  fact  a slight  drop  occurs  at  the  lowest  frequencies.  According  to  the  T1  data, 

r - 1.5  x lo"'  seconds  so  that  T should  rise  around  1 MHz.  It  is  possible  that 
c ^ (17) 

there  are  two  branches  to  the  spectral  density'  \ so  that  the  r , appropriate 

^ /in. 

to  T^s^O-7  seconds,  orthat  extraneous  local  field  effects  are  present'  . Since 
T » T in  a system  which  is  known  to  .be  predominantly  liquid  (cur  « 1)  It  is 

1 2 c (19) 

dangerous  to  infer  anything  from  the  ratio  T./T  In  those  adsorption  systems 

* « (16) 

Fig.  8 shows  the  temperature  dependence  of  Tg  at  full-pores  . Plots 

of  InT  vs.  1/T  are  non-linear  so  that  the  data  cannot  be  fitted  with  a single 
2 

activation  energy.  In  fact  tho  behaviour  is  reminiscent  of  the  behaviour  of  both 

solid  3He*20)  and  a singlo  3He  monolayer  on  Vycor*t2),  with  a 'plateau’  region  occurring 

at  low  temperatures  as  a result  of  exchange  tunnelling.  Various  models  for  Tg 

and  T arc  discussed  in  sections  3 and  4. 

1 (21) 

Experiments  on  the  transport  coefficients  (thermal  conductivity  , mass 

flow*22)  and  spin  diffusion  coefficient*23*)  have  been  carried  out  to  investigate  the 

modification  of  Landau-type  properties  in  the  collisionless  regime.  The  moan- 

free-path  X = .3/T2  (in  A)*24*  and  therefore  becomes  comparable  with  the  pore- 

diameter  d £60  A)  at  approximately  0.1  K.  Rice*25*  has  derived  expressions  for 

the  transport  coefficients  in  the  collisionless  regime,  allowing  for  both  specular  and 

non-specular  boundary  conditions.  For  the  spin-diffusion  coefficient  D,  the  limiting 

expression  (X  »d)  is 
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where  Z is  the  usual  Landau  parameter,  v (he  specular  reflection  coefficient,  V 
o r 

the  Fermi  velocity,  and  d the  pore  diameter  (assuming  perfect  cylindrical  channels). 

All  of  the  transport  experiments  are  in  qualitative  agreement  with  thlB  suggesting 

that  Landau  theory  is  valid  in  the  A.  »d  regime  To  get  numerical  agreement, 

the  data  have  to  be  scaled  by  a factor  F = (X,  where  e is  the  porosity  and  X the 

(26 1 

tortuosity  factor  for  Vycor  . There  is  some  confusion  as  to  the  correct  definition 
of  F'-  and  the  values  used  in  the  three  transport  experiments  differ  considerably'  , 

(29\ 

although  there  may  be  fairly  obvious  experimental  reasons  for  these  discrepancies'  . 
In  the  spin  diffusion  experiment'  ’ F is  determined  by  scaling  the  results  (Fig.  7)  to 
the  high  temperature  region  (~0.2K)  where  A <K  d and  the  particles  diffuse  withir. 
a pore  at  the  ordinary  bulk  diffusion  rate.  F is  taken  to  be  0.  09  + 80  %,  the 
uncertainty  arising  from  the  low  accuracy  of  the  results  around  0.2  K.  Notwith- 
standing some  possible  adjustment  from  the  scaling  factor  there  is  still  discrepancy 
between  data  end  the  predictions  of  Rice  theory  (Fig.  7).  The  shape  of  the  curve  in  the 
A ~d  regime  fe  very  different  from  the  lie  of  the  data,  and  the  experimental  limiting 
value  is  ~half  the  value  given  by  Eq.  (2)  with  v = 0 (non-aero  values  of  v increase 
the  discrpancy).  Rolt  has  Indicated^)  that,  to  explain  certain  features  of  the  data, 
the  adsorbed  layers  should  be  treated  as  an  integral  part  of  the  film  as  a whole.  He 
suggests  that  the  discrepancy  between  data  and  Rice  theory  may  be  related  to  the 
presence  of  the  slow  moving  layers.  It  is  the  purpose  of  the  final  section  of  these  notes 
to  derive  a model  which  takes  account  of  effects  due  to  the  adsorbed  layers  with 
particular  reference  to  the  spin  diffusion  experiment.  References  are  also  made 
to  the  application  of  this  new  model  to  other  transport  experiments. 

2.  THE  NUCLEAR  SUSCEPTIBILITY 


The  ’Two-Fluid1  Model 

For  temperatures  below  0. 1 K the  quasiparticle  coherence  length  (mean- 

3 

free-path)  becomes  comparable  with  the  pore  diameter  allowing  the  He  system 
to  become  increasingly  dominated  by  the  presence  of  the  wall.  Since  (by  definition) 
the  properties  of  a quasi-particle  system  are  determined  by  lnter-quasl-particlo 
collisions,  and  because,  at  very  low  temperatures,  all  lnter-quasi-partlcle 
collisions  occur  at  the  wall,  where  (ignoring  interactions  with  the  first  layer) 
the  quasiparticle  density  is  that  appropriate  to  bulk  liquid  at  32  atm.  the  system 
might  be  expected  to  take  on  'ha  characteristics  of  bulk  helium  at  that  high  density. 
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Justification  for  treating  the  second  layer  region  as  a liquid  stems  from  the 
fact  that  (to  explain  the  spin  diffusion  data)  it  appears  v - 0.  The  relationship 
between  the  second  layer  surface  density  a and  the  sticking  time  T^for  jumps 
Into  the  central  liquid  region  is 

a = (1  - v)  ~ NVp  (3) 

where  N is  the  number  density.  It  is  known  that  a ~ 6.5  x lO^/cm^  ^ \ giving 

T,  ~ 3 x 10  ^ seconds  which  Is  comparable  with  the  'time'  for  exchange  inter- 

actions  to  occur  In  liquid  ' He'  , but  much  shorter  than  the  values  expected  for  a 
-9 

solid  (~  10  seconds,  with  a Neel  temperature  of  36  mdeg) 

These  Ideas  have  led  to  a simple  phenomenological  model  which  we  refer  to 
as  the  ’two-fluid'  model.  It  says  the  following.  Particles  which  collide  in  the 
central  liquid  region  are  ordinary  bulk  liquid  particles  and  they  have  the  properties 
of  bulk  liquid  3IIe  atoms  at  S.  V.  P.  Particles  which  collide  (only)  in  the  second 

3 

layer  region  are  'second  layer'  atoms  and  have  the  properties  of  bulk  liquid  He 
at  32  atm. 

Above  0. 4 K,  X ~ the  interatomic  spacing  so  that  the  'second  layer  type 
particles'  are  confined  to  the  second  layer  region.  In  this  regime  a statistical 
layer  model  is  appropriate  (Eq.  1).  When  X » d the  only  effective  scattering 
collisions  occur  at  the  wall,  so  that  the  'second  layer'  Is  extended  over  the  whole 
pore.  At  Intermediate  temperatures  the  central  region  consists  of  two  types  of 
particles  namely  a fraction  t c which  have  not  collided  since  last  leaving  the  wall 
(to  which  we  ascribe  the  properties  of  layer  2)  and  a fraction  (1  - f)cL  which 
have  collided  and  are  therefore  bulk  particles.  The  system  then  behaves  as  if 
there  were  two  non-interacting  Fermi  systems,  with  some  particles  continually 
traversing  the  pore  and  others  colliding  in  the  bulk  region. 

To  determine  f we  consider  the  following  detailed  balancing  argument.  The 
dominant  Interaction  In  both  bulk  liquid  and  second  layer  regions  is  ferromagnetic 
exchange,  so  th£  an f spin  can  only  leave  the  second  layer!  an  f spin  arrives 
(similarly  for  4 spins).  At  intermediate  temperatures  an  (approximate)  fraction 
f = e particles  leaving  tho  second  layer  at  any  instant  will  traverse  the  pore 
without  collision.  These  particles  have  the  net  polarization  appropriate  to  the 
second  layer  and  cause  an  equal  fraction  with  the  same  polarization  to  leave.  The 
fraction  (1  - f)  with  polarization  appropriate  to  bulk  liquid  (S.V.P. ) Induce  (1  - f) 
similarly  polarized  partlcels  to  leave.  A dynamic  balance  Is  therefore  maintained 
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so  that  at  any  instant  tho  central  region  has  a susceptibility  contribution 

% + “ ' f'*L 

Tho  second  layer  itself  maintains  its  susceptibility  \ y.  Tho  total  film  susceptibility 
is  then 


:2X2  * fCL*2  * (1-[ICIXL 


Fig.  i shows  l iie  result  of  this  analysis. 

(r  reality  ■>'  course  a single  particle  will  eventually  collide  and  objections 
may  be  raised  at  the  suggestion  that  two  non -interacting  systems  are  present. 

However  tho  particles  are  of  course  Indistinguishable!  and  it  fa  just  because  of  tho 
spin-cons<*rving  constraint  of  the  ferromagnetic  interactions  (which  is  one  of  the 
consequences  of  indistlnguishabliUy)  that  it  becomes  immaterial  fas  far  as  tho  net 
polarization  is  concerned)  whether  one  imagines  a particle  traversing  the  pore 
for  ever  or  whether  after  several  traverses  it  collides,  thereby  becoming  a bulk 
atom,  provided  ai  l'io  same  time  another  'second  layer'  particle  is  erected  else- 
where in  the  film.  The  essential  feature  of  the  model  is  that  there  is  a difference  In 
the  result  of  collisions  occuring  in  the  high  density  second  layer  compared  with.  thos>3 
occurring  in  the  bulk  region,  ".’he  known  magnetic  properties  of  bulk  Ho  at  various 
densities'  ' would  tend  to  support  this  contention. 

It  -should  be  clear  that  it  Is  not  suggested  in  the  model  that  all  particles 
leaving  the  second  layer  have  the  characteristics  (i.e.  the  net  polarization)  of  that 
layer.  In  fact  it  foliowis  from  the  detailed  balancing  argument  that  the  fraction 
of  such  particles  Is  i.  The  parameter  f represents  therefore  a Sort  of  quantum 
mechanic rd  accom modal  Ion  coefficient. 


Object  loirs  to  the  Model 

A number  of  objections  have  been  raised  to  the  model''1*1' 

(i;  Given  that  the  dynamic  balance  is  preserved  in  the  model  by  the  ferro- 
magnetic sp'n  conserving  Interaction,  what  turns  up  the  first  spin?  A similar 

3 

problem  arises  In  both  ferromagnetism  and  in  superfluid  He  where  the  Hamiltonian 
has  the  properties  that  if  two  spire  can  pair  then  a lot  more  mnv  as  well*3^. 


The  dctidJa  of  :he  prdrlnp  process  aro  however  related  to  the  existence  of  non-spin 

>nnn  (.m.liif.  Intm.nntlnnd  utn.ll  . . n V hlAn.  flnln n.tMa.i  pirlit>.  nl«M 


conserving  interactions  such  as  X, 


. Non-spln-consorving  terms  are  also 


important  In  any  situation  of  disequilibrium  (following  a DO  pulse  In  a T. 
experiment  would  be  ono  example).  Suppose  In  /Osv  experiment  In  bulk  He  one 
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were  to  cool  suddenly  from  a relatively  high  temperature  where  \ »d,  then 
the  approach  to  equilibrium  would  undoubtedly  depend  on  the  terms  of 
but  the  equilibrium  susecptiblllty  would  be  related  only  to  spin-independent 
terms  such  as  the  interatomic  Van  der  Waal’s  force.  It  should  be  no  surprise 
that  switching  on  a Van  der  Waal's  Interaction  at  the  wall  has  a similar  effect 
on  the  susceptibility 

(ii)  There  is  a feeling  that  one  is  getting  somethingfer  nothing  in  the  model'. 

One  must  remember  that  the  situation  is  rather  different  from  those  high  temperature 
experiments  where  accommodation  effects  are  present.  There  is  no  suggestion 
that  adjacent  'layers'  are  in  thermal  dis-equlllbrium.  The  reason  for  the  drop 
below  the  c lassie ;d  values  of  mean  susceptibility  etc.  is  simply  the  statistical 
consequence  of  indistinguishabillty.  One  might  say  that  the  presence  of  the  wall 
has  the  effect  of  labelling  those  particles  which  collide  with  it,  thereby  Introducing 
an  element  of  distinguishable  with  consequent  removal  of  degeneracy.  (The  two 
fluids  would  then  consist  of  particles  labelled  "wall"  and  particles  labelled  "bulk", 
with  each  collision  rubbing  off  the  old  label  and  creating  a new  one  appropriate  to 
the  region  of  collision'.) 

The  Range  of  Effect 

It  is  very  clear  that  this  model  relates  the  susceptibility  enhancement 

to  a coherence  length  (or  mean  free  path),  so  that  little  enhancement  would  be 

expected  in  a system  large  compared  with  X.  The  question  arises  as  to  whether  some 

(31) 

other  cause  of  size  limitation  might  be  present  . The  susceptibility  enhancement 

requires  a net  increase  in  the  Zeeman  energy  of  the  central  liquid  region.  It  seems 

likely  that  the  film  itself  (probably  the  second  layer)  should  be  required  to  provide 

this  energy  and  the  necessary  entropy.  Assuming  that  the  specific  heat  of  the 

second  alyer  is  equal  to  that  of  bulk  liquid  at  30  atm. , we  find'  ' c/N^kT  ~ 4 

where  N Is  the  number  of  surface  (i.  e.  second  layer)  atoms.  The  magnetic 
s 

enhancement  would  then  be  expected  to  cease  when  the  energy  required  to  turn 
up  N bulk  spins  represented  a significant  imposition  on  the  thermal  energy 
stored  in  the  second  layer.  In  its  attempt  to  polarize  the  entire  liquid  region,  the  second 
layer  is,  according  to  this  argument,  Itself  rendered  more  degenerate  and 
presumably  will  ultimately  become  indistinguishable  from  the  bulk  liquid.  Equating 
energies,  we  find 

‘ - Jlv)" l <t^f» 

F F 


(5) 
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where  T ^ and  T are,  respectively,  the  degeneracy  temperatures  of  bulk 
r r 

liquid  3He  at  30  atm.  andS.V.P.^  K We  then  have 


VNs 


(6) 


This  argument  predicts  the  following:- 

(1)  At  1.6  MHz  and  50  mdog,  N /N  ^ 2x  10*  so  that  one  would  expect 

JL  o 

enhancement  in  Vycor. 

(ii)  A reduction  in  enhancement  as  T ->  0 since  entropies  are  dying  out,  this 

effect  compensating  the  rapid  increase  in  X. 

(34) 

(iii)  A greater  enhancement  at  lower  fields'  ' . 

(iv)  No  appreciable  enhancement  for  the  specific  heat  of  the  liquid  region  (which 
would  be  open  to  experimental  study). 


Subsidiary  Comments 

(1)  There  is  an  obvious  analogy  between  the  X »d  regime  in  Vycor  and  the 
electron  system  in  metals,  with  the  substrate  potential  playing  the  part  of  the 
lattice  interaction.  Just  as  a lattice-mediated  electron-electron  coupling  can 
give  rise  to  superconductivity,  it  should  be  possible  for  a substrate-mediated 
interaction  to  produce  (or  enhance  the  possibility  of)  superfluidity  in  3He^35^. 

One  can  easily  Imagine  how  a metastable  state  might  be  nucleated  at  the  walls 
where  the  pressure  is  high.  The  stability  of  this  state  would  presumable  be 
related  to  the  ability  of  the  film  to  polarize  the  liquid.  Bearing  in  mind  the  way 

3 

in  which  the  degree  of  polarization  of  He  depends  on  pressure,  temperature 
and  field*  the  application  of  equation  5 would  predict  a phase  diagram  qualit- 

3 

atlvely  similar  to  that  seen  in  superfluid  He  at  2 mK’.  One  unusual  feature  has 

laready  been  seen  in  superfluid  He  confined  to  ~12  p platinum  powder'  . 

It  is  found  that  the  susceptibility  drops  as  one  cools  into  the  A-phase,  whereas 

In  bulk  a slight  rise  Is  seen*36\  While  there  may  be  other  explanations  of  this 
(37) 

phenomena  ' , it  could  be  that  the  normal  susceptibility  is  enhanced  by  the  effects 
suggested  in  this  model,  but  that  in  the  superfluid  state  the  mechanism  breaks 
down  allowing  the  susceptibility  to  relax  back  to  its  true  bulk  value.  Since  the 
experiments  were  done  at  fairly  low  resonant  field  (50  gauss)  some  enhancement 
might  be  still  expected  despite  the  larger  pore-size  (equation  5). 

(2)  Similar  polarization  effects  are  seen  in  other  systems  where  long  range 
correlations  are  present.  One  such  is  the  proximity  effect  of  superconductivity' 


J 
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where  a material  below  its  own  transition  temperature  can  Induce  superconductivity 
in  a sample  of  normal  metal  simply  by  contact.  Examples  of  ferromagnetic 

materials  inducing  ferromagnetism  in  materials  above  their  owi  Curie  temperature 

139)  3 

have  also  been  seen  , and  in  this  context,  one  must  remember  that  bulk  Ho 

to  very  close  to  becoming  a forromagnet'  . 

(3)  There  is  a number  of  obvious  simplifications  in  the  proposed  model  which 

I 

may  contribute  to  the  fact  that  agreement  between  model  and  data  is  not  perfect  (Fig.  1) 

-X/X 

(1)  The  parameter  f should  really  be  derived  by  Integrating  e over  all 

path  lengths  X allowing  for  a cosine  distribution  of  desorption.  Calculations 

have  shown  that  by  ignoring  this  we  have  underestimated  y^  by  ~2%.  In  view 

of  the  uncertainty  in  the  exact  nature  of  the  Vycor  geometry  with  the  possible  existence 

of  pores  much  smaller  than  the  average  (which  would  be  strongly  weighted  in  the  f factor) 

we  have  not  attempted  too  sophisticated  an  analysis. 

(li)  We  have  Ignored  assimilations  of  bulk  particles  directly  into  the  first 
layer.  This  seems  justified  in  view  of  the  oxpected  correlation  time  for  inter- 

_7 

layer  motion  (believed  to  £ 10  secs)  but  the  effect  of  such  exchanges  will 
become  significant  at  very  low  temperatures  where  the  monolayer  susceptibility 
will  be  many  orders  of  magnitude  greater  than  the  bulk  liquid.  To  explain  the 
discrepancy  between  model  and  data  (at  40  mK)  in  terms  of  this  effect  alone  requires 
approximately  1 in  20  incident  particles  to  penetrate  into  the  first  layer.  Since 
Tj_  for  the  socond  layer  is  ~ 3 x 10  ^ seconds  this  would  give  t.  for  the  first 
layer  the  unlikely  value  of  6 x 10  seconds,  much  shorter  than  the  value  (10  secs) 
believed  to  bo  correct  (Section  1). 

(ill)  A slight  modification  of  the  model  has  been  examined,  allowing  for  the  fact 
that  in  the  first  collision  after  leaving  the  wall,  residual  wall  effects  may  still  be 
felt  as  the  quasiparticlc  wave  packet  is  still  coherent  with  the  wall  interaction 
Only  after  the  first  collision  is  this  coherence  lost.  The  'second  layer'  type 
particles  would  then  bo  identified  with  those  particles  which  have  had  not  more 
than  one  collision.  Again  an  increase  in  y^  is  expected. 

(iv)  We  have  ignored  the  extended  density  profile  which  was  found  necessary  to 
explain  the  high  temperature  susceptibility  data^\  The  'third'  layer  should  be 
accredited  with  a susceptibility  slightly  greater  than  that  of  bulk  liquid  at  the 
S.  V.  P.  Calculations  show  that  we  have  again  underestimated  by  ~3%. 

All  of  these  simplifications  tend  to  underestimate  y^  and  therefore  it 
should  be  fairly  easy  to  reconcile  the  residual  discrepancy. 
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An  alternative  picture 

In  the  previous  model  the  first  layer  was  essentially  ignored,  although,  as 
already  mentioned,  it  could  be  easily  indided  (and  perhaps  should  be  so  included 
at  very  low  temporatures).  However,  since  all  of  the  spins  are  exchange-coupled, 
it  is  perhaps  better  to  treat  the  entire  film  as  a unique  three-dimensional  entity 
:ind  to  attempt  to  predict  the  wavefunctions  for  this  system. 

The  standing  wave  eigenfunctions  for  a 'particle  in  a box'  are  well  known. 

If  one  now  switches  on  a highly  localized  but  large  binding  interaction  at  the  walls, 
the  immediate  effect  is  to  me  l ate  particles  (i.e.  increase  certain  single  particle 
wavefunctions)  near  the  wall.  However  having  allowed  far  hard-core  repulsion 
(in  a self-consistent  field  calculation),  the  attractive  Interaction  seen  by  incident 
particles  as  they  approach  the  wall  becomes  a repulsion,  and  the  average 
potential  seen  by  a particle  looks  rather  like  and  S.H.O.  potential**0^.  In  such 
a system  the  thermally  active  high  quantum  number  particles  near  the  top  of  the 
Fermi  distribution  will  have  wavefunctions  which  are  macroscopically  large  near 
the  walls,  so  that  once  again,  the  behaviour  of  the  relevant  particles  will  be 
dominated  by  the  wall  interaction. 

In  a sense  this  picture  is  predicting  a 'condensation'  of  the  active 

particles  into  a region  near  the  wall, while  simultaneously  conserving  overall  mass  by 

allowing  the  wavefunctions  of  those  particles  lower  down  the  momentum  distribution 

to  have  wavefunctions  peaking  in  the  central  region.  The  relationship  betweeen 

this  typo  of  picture,  and  the  idea  that  the  second  layer  is  solidifying  (or  that  its 

(7) 

correlation  time  is  increasing)  as  the  tomperaturo  Is  reduced  , is  self- 
evident.  Since  it  seems  unlikely  that  one  will  be  ever  able  to  treat  a system  of 
hard-core  particles  and  then  make  a quasi-particle  transformation^,  the  ad- 
hoc  models  probably  deserve  closer  attention. 

This  typo  of  picture  is  similar  to  the  mean-free-patb  (two-fluid)  model 
in  that  enhancement  effects  are  related  to  the  behaviour  of  the  quasiparticle 
wavefunctions  in  the  first  and  second  layer  regions.  However,  the  big  difference 
is  in  the  range  of  effect  predicted.  In  a larger  system  the  average  potential 
seen  will  be  affected  much  less  by  the  wall  Interaction  and  It  Is  unlikely  that  a 
build-up  of  wavefunction  could  occur  in  the  same  way.  Measurements  on  adsorption 
systems  with  larger  pore-size  would  therefore  allow  one  lb  discriminate  hntwoon 
these  alternative  explanations. 
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3.  SPIN-SPIN  EQUILIBRIUM  IN  THE  SATURATED  FILM 

To  simplify  a discussion  of  tho  dephasing  mechanism  in  the  saturated 
film  we  begin  by  using  the  layer  model.  This  is  obviously  an  abstraction  at 
low  temperatures  where  the  film  behaves  as  a three-dimensional  entity. 

The  Dephaslng  Zone 

T^  for  the  full*-pores  film  is  very  much  shorter  than  the  T 2 of  bulk 

liquid  and  comparable  (slightly  longer)  than  the  monolayer  value^^. 

The  correlation  time  Tu  for  motion  within  a single  monolayer  is  believed 
-7  * 3 

to  be  10  seconds.  A solid-like  treatment  of  a second  layer  of  He  (adsorbed 

4 -9 

on  He)  has  suggested  that  tc  is  of  order  10  seconds.  Allowing  for  a return 

towards  liquid-like  behaviour  when  that  layer  is  depleted  on  subsequent 

adsorption  , it  is  concluded  that  for  jUl  layers  above  the  first  the  correlation 
-9 

time  is  < 10  seconds  and  therefore  that  T„  » several  milliseconds.  We 

mt 

identify  , therefore,  the  first  layer  as  the  dephasing  zone  and  assume  that  all 
spins  dophase  by  Interaction  with  (or  assimilation  into)  that  layer. 

T2  will  then  be  related  to  T,^2' the  spin-spin  oquilibrium  appropriate 
to  the  monolayer,  and  to  the  relative  amount  of  time  r.  the  spins  spend  there. 
Since,  in  a time  t,  particles  will  make  t/cr/c^)  visits  to  the  zone  and  spent  c^t 
seconds  relaxing,  tho  time  scale  is  effectively  expanded  by  the  ratio  c , and 
this  would  be  consistent  with  an  approximate  factor  of  two  in  the  ratio  T /T^  ’ 

(figs.  4 and  5).  However,  the  exact  relationship  depends  on  a number  of  factors. 
Four  situations  can  be  envisaged**1*. 

The  Four  Regimes 

(i)  If  the  time,  t , for  particles  to  diffuse  to  the  wall  Is  » than  either  T *2* 

w ^ -9 

or  t,  . then  t is  tho  controlling  factor.  In  the  Vycor  system  this  time  Is  ~ 10 

X u) 

seconds  so  thrt  this  regime  does  not  apply. 

(11)  If  the  specular  reflection  coefficient  for  assimilation  Into  the  first 

layer  Is  high,  so  that  particles  make  many  traversals  between  second  layer  and  bulk 

before  assimilation  occurs,  then  t,  may  become  the  controlling  factor.  Provided 
^ (l)  . 

Tl  * T2  we  expect 

T2  = r/Cj  (7) 
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(ill)  The  regime  t ^r,  T ' has  two  Bubreglons. 

uJ  * 2 

(a)  - Ti; ' w^ore  th  18  the  correlation  time  for  motion  within 

the  first  layer. 

In  this  situation  a particle  visits  a layer  for  time  r and  sees  a local  field 
i 

M 3 which  is  virtually  unchanging  over  the  length  of  its  stay.  During  this  visit 

2 h 

it  dephases  by  an  amount  A#  - tM  3.  The  particle  then  leaves  the  zone  for 

2 l 

a time  (1  - c )/tc  , and  then  re-enters  to  find  a local  field  of  magnitude  M_z 

i 1 Ct 

but  probably  of  opposite  polarity,  so  that  it  Cephases'  in  the  opposite  direction. 
Avoragod  over  many  cycle  times  ( r/c ) the  particle  makes  tc Jr  visits  and 
dephases  by  A0(t)  where 


- 


L — J A* 


18) 


The  parameter  f is  a numerical  factor  of  order  unity  which  can  be  related  to 
the  spectral  density  by  the  equation 


‘(wt)  = [1  + ^rJ(wr)  + |j(2wf)] 


(9) 


2 2 

where  typically  j(wr)  = 1/(1  + w r ).  By  definition,  T 2 is  the  time  for 
which  A0  = 1,  giving 


T - i_  J_  _i_  _i_  T (®) 

2 cl  f2  M2T  C1  2 * 

The  function  f is  of  course  chosen  to  reflect  the  frequency  dependence  of  T ^ 
and  will  display  the  usual  10/3  effect.  This  analysis  shows  that  if  interactions 
near  the  monolayer  are  isotropic  (t^  ~ r^)  then  is  longer  than  T ^ t*ie 
relative  amount  of  time  spent  in  the  dephaslng  zone. 

(b)  Tg  '»  t^»  Tg  . During  a single  visit,  a particle  now 
experiences  (T^/Tg)  Jumps  and  debases  by  an  amount 


/TJL  n2  £ 

< A 0 > = f(«T„)  (— ) M.  i 


2 T« 


(10) 


(11) 


having  visited  a whole  spectrum  of  sites.  After  leaving  for  time  (1  - 
it  re-enters  the  zone  and  again  dephases  by  < A0  > . In  time  t » it  will  have 
made  (tc^/r^  ) visits  and  therefore  dephased  by 

— * ««-V 


<A0> 


(12) 
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giving  the  following  expression  for  T , 

M 


j_  _i_  /Ik.  , i (T  wj 

°i  "“V  ' T»  J M.i  ci  2 


8 2 (4) 

To  fit  the  data  with  a second  moment  M = 4 x 10  (rad/sec)  ' , r.  must  be 
-5  -4  2 

approximately  10  to  10  seconds  which  Is  consistent  with  the  necessary 

-3 

requirement  ~ 10  seconds. 

The  Temperature  and  Field  Dependence  of 

We  now  examine  these  models  in  relation  to  the  variation  of  T with 

temperature  and  field.  Changes  might  be  expected  in  , c^,  and  T^ 

If  is  the  controlling  factor  (model  £1))  then  r must  be  ~ several 

milliseconds,  and  must  Increase  at  higher  temperatures.  It  is  difficult  to  Imagine 

a situation  in  which  could  increase  with  temperature,  even  if  both  quantum 

tunnelling  and  thermal  activated  motion  are  present.  If  one  can  justifiably  take 

c.  to  be  proportional  to  r.  (as  in  equation  3)  then  the  rather  slow  field  dependence 

(at  0.1  K)  could  be  understood.  Since  r ~ T ' ' the  temperature  dependence  might 

(z)  2 

be  complicated  by  the  rise  in  T ' ' with  temperature,  so  that  at  higher  temperatures 
regime  (ill)  becomes  appropriate. 


If  situation  (lll)(a)  is  operative  the  increase  of  T^  with  temperature  could 
be  attributed  directly  to  the  variation  of  T The  field  dependence  is  however 

difficult  to  explain.  T^  for  a monolayer  increases  monotonic  ally  with  frequency, 
through  the  f factor,  whereas  T for  the  saturated  film  shows  a very  slow  drop 
with  increasing  field.  One  can  attribute  the  difference  only  to  an  Increase  in  ^ 
with  increasing  field.  At  higher  fields,  therefore,  one  muBt  move  into  regime 
(iii)(b)  where  the  term  T„^  enters  as  a square-root.  The  monotonlcally  increasing 
behaviour  cannot  be  rendered  decreasing  simply  square -rooting.  The  difference 
would  be  ascribed  to  changes  in  the  ratio  T^/c. . The  inferences  one  can  draw 
from  this  now  depend  on  whether  c^  « or  whether  c^  and  r^can  vary  independently. 
The  former  would  require  an  increase  in  T^(and  therefore  c^),  the  latter  a drop 

in  t (with  c.  constant).  Arguments  can  bo  put  forward  for  both  possibilities. 

• 1 3 uv-2) 

On  the  basis  of  the  field  dependence  of  the  bulk  He  melting  curve*^,  a length- 


ening of  might  be  expected.  Allowing  for  the  polarizing  effects  discussed 
in  equations  5 and  6,  the  opposite  might  be  expected. 

The  temperature  dependence  of  T on  fills  model  would  also  be  consistent 
with  r.  getting  longer  as  the  sample  is  cooled.  This  may  be  evidence  of  a 
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progressive  condensation  (solidification)  of  the  second  layer  onto  the  first  layer^. 

Alternatively  it  may  reflect  the  increasing  dominance  of  the  monolayer  over  the 

entire  system.  Since  in  any  magnetic  experiment  the  three  components  of  the  film 

are  weighted  according  to  their  contributions  to  the  susceptibility  (rather  than 

to  the  relative  number  of  spins  in  each  region) 

it  may  be  more  correct  to  replace  c^  by  c^/)^  whore  ^ Is  the  (normalise^ 

susceptibility  of  the  whole  film  (Fig.  1).  Expressions  10  and  13  then  become 

(assuming  c “ T );- 
* 


- c7 

<TX  ~ V 

(10’) 

and 

T2 

* V 

(13*) 

At  high  tempe  rat  ires,  of  course,  the  respective  equations  are  identical.  The 
data  cannot  be  fitted  exactly  to  equation  10*  over  the  entire  temperature  range 

although  the  order  of  magnitude  variation  is  correct.  However,  up  to  0.15  K 

1 

the  ratio  T^/)^  is  constant  to  within  2 %,  suggesting  that  equation  13'  might 
be  correct  with  a temperature  independent  T ^ = 1.5  msecs.  The  faster 
variation  of  T above  0. 15  K is  explicable  on  two  grounds.  Either,  as  r, 
shortens,  one  moves  in  regime  (lii)(a)  described  by  equation  10',  or  Tg^ 
ltsolf  starts  to  increase  (Section  I).  Measurements  at  lowet  temperatures 
would  be  useful  in  exploring  the  ratio  Tg/)^  . 

All  of  the  models  (il),  (iii)(a),  (iii)(b)  are  plausible.  To  discriminate 
between  them  more  accurate  measurements  of  T 2 should  be  made  over  larger 
frequency  and  temperature  range  both  for  full-pores  and  a single  monolayer. 

The  detailed  relationship  between  T and  T ^ might  then  be  established.  It 

2 2 3 

would  also  be  instructive  to  measure  T for  a single  layer  of  He  adsorbed  on 

4 

He  - on  the  basis  of  these  arguments  T is  expected  to  be  long  (~  1 second). 

If  c ^ is  itself  field-dependent,  the  nuclear  susceptibility  may  bot  be  uniquely  defined. 
Measurements  of  \ against  temperature  over  a range  of  frequency,  scaled  to  the 
high  temperature  region,  would  be  particularly  useful. 

-7  -3 

Whichever  model  is  correct,  it  seems  clear  that  10  < r & 10  seconds. 

X 

It  also  seems  likely  that  -r  must  change  with  field  and  temperatire  to  explain  tte  T 

it 

behaviour. 


In  this  analysis  we  have  assumed  that  effects  due  to  extraneous  fields  (im- 
purities etc. ) were  either  absent,  or  present  to  the  same  extent  in  both  full-pores  ami 
monolayer  experiments,  ff  the  latter  is  the  case,  it  remains  unlikely  that  the  very 
different  frequency  dependence  of  T could  be  explained  in  terms  of  extraneous  effects 
alone. 

The  Two-Fluid  Model 

In  the  analysis  of  T9  given  above  we  have  assumed  that  the  dephasing  zone 
stays  'intact*  as  more  layers  are  added,  ff  this  is  so,  then  it  makes  sense  to  examine 
the  ratio  T /T,  At  low  temperatures,  however,  it  is  probably  more  correct  to 
treat  the  entire  film  as  a three-dimensional  entity  and  allow  for  the  fact  that  the  emo- 
tional behaviour  in  the  zone  may  be  completely  transformed.  The  temperature  and 
frequency  response  of  T ^ itself  may  therefore  be  very  different  from  that  of  a single 

S 

monolayer. 

A 'two-fluid-model'  type  of  approach  has  been  made,  based  on  the  equation 
-1  m <2> 

T2  = f/T2  + <1_f>/T2  (14) 

and  T92\  which  have  been  treated  as  temperature-independent  parameters,  rep- 
2 2 

resent  the  spin-relaxation  times  of  the  two  components.  Since  the  'second  layer'  type 
particles  are  the  ones  which  are  dominated  by  the  presence  of  the  adsorbed  layers  we 

expect  T„^  to  be  less  than  T9?'. 

« ^ (2\ 

This  analysis  has  shown  that  (1)  for  no  values  of  T9  ' and  T'  can  the  data 

be  fitted  to  equation  (14)  with  f = c K and  d = 60  (ii)  an  almost  perfect  fit  is  ob- 
tained with  d = 1 0 % and  T^  * 1.6  msecs,  and  T^  = 2.9  , secs. 

The  interpretation  of  (ii)  Would  presumably  be  that  between  the  antlforromagnetic 
(dephasing)  region  and  the  ferromagnetic  core,  an  Intermediate  zone  of  width  10^ 
exists,  with  the  property  that  only  collisions  within  this  region  can  cause  spins  to  lose 
coherence  with  the  dephasing  centre.  Since,  with  three  variable  parameters,  a good 
fit  to  the  data  is  bound  to  occur,  it  will  be  necessary  to  make  more  precise  measure- 
ments in  large  geometries  to  evaluate  this  idea. 

4.  SPIN- LATTICE  RELAXATION  IN  3He  FILMS 

The  significant  feature  of  the  T1  data  is  that  T1  « f2,  characteristic  of 
three-dimensional  fluids,  while  the  values  of  Mg  and  T^  are  appropriate  to  the  mono- 
layer density  (section  1).  Since  the  system  is  behaving  in  this  three-dimensional 
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way  It  is  probably  wrong  to  attempt  to  relate  to  the  T ^ of  a single  monolayer. 

However,  bearing  in  mind  that  the  second  and  higher  layers  have  very  short  correl- 

3 

ation  times  so  that  T should  be  long  and  frequency-independent  (T1  for  a layer  He 

on  He  was  '■Jsecond  at  6 MHz)  it  is  clear  that  over  the  range  of  frequency  on  Fig  3 , 

(43) 

the  monolayer  region  is  the  relaxation  zone v ’ . 

The  nature  of  the  assimilation  process  is  not  understood  in  detail  but  must 

relate  to  the  interlayer  exchange  motion  (rather  like  the  diffusion  of  exchange  energy 
3 (44) 

in  solid  He'  ').  The  following  points  seem  relevant. 


1.  Fluctuation  Effects 

Immediately  after  a 90*  pulse  the  spin  populations  in  the  central  region  are 

equalised.  If  an  incident  4 particle  penetrates  the  antiferromagnetic  relaxation  zone, 

then  an  t spin  will  be  injected  into  the  central  liquid  region.  Effectively  the  liquid 

has  taken  one  step  towards  returning  to  its  equilibrium  polarization  and  the  4 spin, 

while  in  the  relaxation  zone,  stands  some  probability  of  being  'flipped-up'  itself.  One 

must  remember  that  the  newly  released  spins  will  eventually  return  to  the  zone.  If 

this  occurs  in  a region  which  has  Just  experienced  a net  recovery  ( 4 to  t ) then  an 

antiferromagnetic  interaction  with  an  incident  particle  will  reverse  the  net  polarity 

of  this  region.  We  therefore  envisage  each  region  of  the  zone  fluctuating  rapidly  in 

net.  polarity,  with  the  average  rate  of  advance  from  ^ to  f (caused  by  the  relaxation 

l/2 

process)  proportional  to  (the  number  of  visits)  . Assuming  the  probability  of  re- 
covery in  a given  visit  (fluctuation)  to  be  history-independent,  the  relation  for 
becomes 


i 


(1  - 1/e)  « 0. 6 


elvln*  T,  - (tJV/tj.  H5> 

This  equation  would  explain  therefore  the  relationship  between  the  frequency  depen- 
dencies of  the  full-pores  film  and  the  monolayer.  However,  there  are  two  features 

which  are  unsatisfactory.  First,  to  give  numerical  agreement  we  must  have 

-6  -3 

~ Tj  ~ tens  of  milliseconds,  whereas  the  data  indicate  that  10  £ r jg 10 

seconds.  Second,  equation  (15)  suggests  that  T1  for  the  saturated  film  should  be  linear 
through  the  origin. 


- 22  - 


2.  Uncorrelated  Interlayer  Motion 

Since  the  excahnge  Interaction  conserves  spin  it  cannot  itself  produce  relax- 
ation and  it  is  not  obvious  that  interlayer  exchange  motion  can  ever  produce  a net 
change  in  the  polarization  of  particles  injected  into  the  central  region.  The  relax- 
ation process  may  depend  on  (uncorrelated)  single-particle  quantum-tunnelling  into 
the  zone.  Assuming  visits  to  the  zone  are  statistically  independent,  we  may  write 

T1  = T1Z>  1 Cl  * (16) 

To  explain  the  different  frequency  and  temperature  dependence  In  terms  of  c^  alone 
(ie  with  ^qualtothe  monolayer  value)  the  behaviour  of  c ^ would  have  to  be  very 
different  from  that  expected  from  considerations  of  T0  (Section  3).  For  example 
would  have  to  decrease  with  frequency,  and  to  be  essentially  independent  of  temper- 
ature (Section  1). 


3.  The  Removal  of  Anisotropy 

The  appearance  of  an  intercept  in  the  T^  v.  f^  plot  suggests  that  the  cause  of 

the  unusual  feature  of  the  monolayer  (T^  v.  frequency  passing  through  the  origin)  has 

been  removed  in  the  multilayer  film.  If  that  cause  was  surface  heterogenlty  then, 

since  spin-lattice  relaxation  occurs  as  a result  of  'assimilation'  into  the  monolayer, 

the  range  of  correlation  times  should  still  have  effect.  An  explainatlon  can,  however, 

be  put  forward  on  the  basis  of  the  removal  of  anisotropy  effects. 

In  a single  monolayer  spins  are  constrained  to  lie  and  move  on  surfaces  at 

fixed  angles  to  the  field  axis.  On  some  surfaces  those  components  of  KdiP  c,,Bble 

of  inducing  relaxation  are  identically  zero.  Having  made  many  circumnavigations 

of  a pore  circumference  in  time  ~ T^,  a spin  comes  to  two  conclusions'.  First  It 

sees  2.  mean  value  of  < X ..  > averaged  over  all  surfaces.  Second,  it  'remem- 

dip 

bers1  that,  on  some  surfaces,  it  had  to  wait  a lot  longer  for  3C^  *°  c*,an8e's^8n^" 

lcantly  and  it  therefore  concludes  that  a range  of  r *s  Is  present.  It  is  shown  else- 

c 

where  that  a range  of  re's  can  lead  to  a linear  frequency  dependence  (12,  14).  IT, 

In  the  multilayer  film,  motion  near  the  monolayer  region  Is  sufficiently  Isotropic 
<T  ^ Tjj)  then  this  constraint  will  be  removed  And  the  system  described  by  a single 
correlation  time. 
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4.  Extraneous  Field  Effects 

There  is  no  obvious  way  of  explaining  the  different  frequency  dependencies  of 

T in  the  monolayer  and  multilayer  films  in  terms  of  a field-dependent  'M  ’.  Further- 

* 8 2 ^ 
more  the  value  M ■ 2.  5 x 10  (rad/sec)  derived  from  Fig  3 is  that  expected  on  the 

basis  of  intrinsic  interactions  alone. 

8 2 

One  problem  still  unresolved  is  that  if  M * 2.  5 x 10  (rad/sec)  and 

fj  2 

r = 1.5  x 10  1 seconds,  then  T„  should  drop  for  frequencies  below  JMH  . 
c 2 • z 

The  constancy  of  T (Fig  5)  may  be  due  to  a low  frequency  ’hump'  In  the  spectral 

2 (45) 

density  which  is  affecting  but  not  . 

5.  The  Approach  to  Equilibrium 

One  of  the  most  interesting  problems  concerns  the  shape  of  the  recovery 
in  time  with  particular  reference  to  the  way  the  susceptibility  returns  to  its  equili- 
brium value  in  which  the  susceptibility  profile  is  maintained.  Since  relaxation  is 
occurring  in  the  highly  polarized  monolayer,  the  system  may  tend,  in  the  early  stages 
of  recovery,  towards  the  ’Curie'  susceptibility,  before  approaching  its  degenerate 
equilibrium  valve  with  a slower  time  constant.  Non-exponential  recoveries  of  this 
form,  are  seen^4®\  At  0.4  K the  system  is  ~ 30%  degenerate  in  a regime  wher?  the 
pore-diameter  is  many  orders  of  magnitude  smaller  than  the  mean-free- path  for 
non-spin-conserving  (relaxative)  interactions.  It  follows  that  the  enhanced  polari- 
zation 8<*en  ?t  lower  temperatures  cannot  be  related  to  dynamic  polarization  effects 
alone.  The  'two-fluid'  model  of  course  relates  the  enhancement  to  the  number  of 
spin  consex  ving  interactions,  so  that  when  X » d surface  polarization  effects  may  be 
possible. 

A conclusion 

The  most  self-consistent  explanation  of  the  T^  data  is  as  follows.  The  mul- 
tilayer film  behaves  as  three-dimensional  fluid  with  a lorentzian  type  spectral  density. 

-7 

Anisotropy  effects  seen  in  a single  layer  are  removed.  T = 1. 5 x 10  sees  and 
8 2 

M *a8. 5 x 10  (rad/sec)  . Relaxation  processes  occurring  in  the  monolayer  region 
are  communicated  to  the  rest  of  the  film  via  unspecified  quantum  mechanical  motion. 

Definitive  experiments  are  likely  to  involve  (i)  the  removal  of  the  relaxation 
zone  by  a layer  of  4He  (ii)  separation  of  the  irherent  and  polarizing  properties  of  the 
first  layer  in  larger  geometries  (provided  t^  « Tj). 
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5.  THE  SPIN  DIFFUSION  COEFFICIENT 

The  basic  idea  behind  the  transport  coefficient  experiments  in  the  X >>  d 
regime  has  already  been  discussed.  These  experiments  tend  to  show  that  the  o-;'.sl- 
particle  concept  is  valid  in  restricted  geometries  and  the  specular  reflection  coef- 
ficient v a 0.  The  theories  of  Rice  are  approximately  obeyed  although  the  spin 
diffusion  coefficient  is  smaller  than  the  expected  value  (Fig  7).  We  examine  below 
in  more  detail  the  role  the  adsorbed  layers  might  play  in  the  spin  diffusion  experi- 
ment. These  effects  which  were  completely  Ignored  by  Rice  are  expected  to  be  large 
at  low  temperatures  where  the  monolayer  susceptibility  is  very  much  greater  than 
that  of  the  liquid.  There  are  three  approaches  all  of  which  are  derivatives  of  the 
two-fluid  model.  First  we  allow  for  the  polarizing  effect  of  the  layers.  Next  we  con- 
sider the  delaying  effect  experienced  by  particles  when  they  are  adsorbed.  Finally, 
we  Include  the  effects  of  surface  diffusion. 

The  Two-Fluid  Model 

The  idea  that  the  system  behaves  as  if  two  non-interacting  *fluids ' is  present 
has  led  to  a reasonably  self-consistent  description  of  the  net  polarization  of  the  film 
(Section  2).  It  is  somewhat  more  difficult  to  Justify  this  model  when  applied  to  trans- 
port experiments,  but  the  picture  preserves  the  essential  physical  feature,  namely 
that  those  particles  with  higher  susceptibility  (the  second-layer  ones)  are  also  the  ones 
which  travel  farthest.  We  ignore  motion  in  the  adsorbed  layers,  treating  the  central- 
liquid  region  as  the  only  mobile  part  of  the  film.  This  may  be  conceptually  correct 
because  the  dominant  Interaction  is  ferromagnetic  exchange,  so  that  it  is  Immaterial 
whether  one  imagines  one  incident  quasip.rticle  to  be  adsorbed  and  another  randomly 
desorbed,  or  whether  one  particle  remains  in  the  layer  and  merely  incoherently  scat- 
ters the  incident  particle. 

Since  there  are  two  non-interacting  fluids  (one  in  the  Knudsen  regime,  the  other 
(47) 

in  collioional  regime)  , we  add  the  diffusing  coefficients  of  each  component  (rather 
than  adding  scattering  rates,  which  are  proportional  to  D 1).  The  total  diffusion  coef- 
ficient D will  be  given  by 

D - QDfc  ♦ (1  - c.)PL  (17) 

where  c * f X2  / <*  Xg  ♦ (1  - 0 /L  ) (18) 


f 
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and  f = e (X  being  the  mean  free  path  of  the  total  system  and  therefore  equal  to 

3 

X of  bulk  He).  The  Kwudsen  term  will  be  taken  as  the  limiting  expression  given  by 
Rice,  namely 


DK  - 1/3  (1+V4»Vc 


(19) 


where  dy  is  the  controlling  pore  diameter.  There  is  considerable  evidence  that  a 
factor  of  two  variation  in  pore-size  is  present  in  Vycor  and  we  take  d = 36^, 

/I  ^ 

half  of  the  volumetrically  determined  radius'  . To  relate  to  the  experimental  values, 

the  scaling  f actor  F = cX  has  to  be  used  and  we  take  F C. 09,  the  value  used  in  the 
(23) 

original  work'  . D^,  the  bulk  liquid  diffusion  coefficient,  is  ascribed  to  those 
particles  colliding  in  the  central  region. 

(491 

Fig  7 shows  D calculated  on  this  model  and  quite  good  agreement  is  obtained  . 
In  calculating  we  have  used  1 + Zq/4  = 0. 26,  the  value  appropriate  to  30  atm. 
liquid  (this  may  be  open  to  objection  but  the  value  for  S.  V.  P.  is  not  very  different, 

0. 33  in  fact).  The  specular  reflection  coefficient  was  assumed  to  be  zero.  Values 
of  v greater  than  zero  would  worsen  agreement  between  experiment  and  model. 

3 

Since  number  density  is  conserved  in  the  central  region  we  assume  V = 5. 38  x 10 

F 

cm/sec  . The  value  of  DK  Is  then  17  x l(f5  cm2/sec.  Sizeable  contrib- 
ution apparently  comes  from  the  few  collisions  that  do  occur  in  the  central  region  at 
low  temperatures  and  it  may  be  that  this  relates  to  the  fact  that  a single  collision  can 
scatter  particles  down  a pore  (with  little  subsequent  probability  of  colliding)  much 
more  effectively  than  the  Knudsen  process. 

New  Model  Allowing  for  Delaying  Effects 

The  main  objection  to  the  previous  model  is  that  to  get  agreement,  the  value 
d^  * 36  R was  chosen  (with  however  some  justification  The  fact  that  the  low 

temperature  data  are  smaller  than  expected  may  be  related  to  the  presence  of  the  slow- 
ly  moving  adsorbed  layers'  . In  the  absence  of  a microscopic  theory  (or  the  lik— 
lihood  of  one)  we  present  a phenomenological  model  allowing  for  the  delaying  effect 
of  the  quasiperticle  wavefunction  when  particles  are  adsorbed.  We  shall  assume  that 
the’ life-history’  of  a single  quasi-particle,  as  it  travels  down  a pore,  can  be  replaced 
hy  an  ensemble  average  of  three  non-interacting ’fluids’,  namely  particles  exchanged 
between  the  first  layer:  particles  exchanged  between  the  second  layer:  and  particles 
colliding  in  the  bulk  region.  The  relative  proportions  of  these  fluids’  will  of  course 
be  highly  temperature  dependent.  Because  of  the  uncertainty  in  the  value  of  the 
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scaling  factor  F (section  1)  we  shall  allow  F to  vary  within  its  range  of  experimental 
error  and  choose  the  valve  shich  gives  the  best  agreement  between  model  and  experi- 
ment. 

The  main  problem  in  any  kinetic  model  of  the  film  structure  is  that  of  quan- 
tifying the  accommodation  coefficient,  the  parameter  which  determines  the  fraction 
of  particles  desorbed  from  a layer  with  the  characteristics  of  that  layer.  We  assume 
on  the  basis  of  the  discussion  in  (section  2)  that  this  coefficient  is  f = e”^*.  In 
other  words,  the  proportion  of  particles  leaving  a layer  with  the  characteristics  of 
that  layer  is  equal  to  the  number  of  similar  particles  arriving.  Assimilation  in  the 
sense  of  this  model  is  a ’quantum  assimilation’,  the  reasons  for  the  different  proper- 
ties in  each  layer  being  quantum  mechanical  in  origin. 

The  assumptions  of  the  model 

a)  The  number  of  particles  incident  on  any  surface  N = 1/4  n V 

r 

b)  A fraction  f Nof  these  come  from  either  1st  or  2nd  layer  without  collision 

-d/A 

where  f = e 

c)  gfN  particles  will  penetrate  into  the  1st  layer,  g is  treated  as  a variable 
parameter  but  is  likely  to  be  small  for  reasons  given  in  d). 

d)  the  surface  density  a for  any  layer  is  related  to  the  sticking  time  r by 

f--  N r ' for  the  second  layer 

~ 5 x 10  ^ seconds  ; 

-10 

~ 10  seconds. 

tc  (for  motion  within  the  first  layer)  is 

„ 10  seconds,  ff  the  layer  model  is  to  be  meaningful  then  j.  (the  Inequality 

1 c 

allowing  for  an  interchange  between  1st  and  2nd  layer  without  an  incident  quasiparticle 
from  the  liquid  region).  Hence  g <<  1. 

e)  the  mean  time  T for  a particle  to  pass  a distance  u down  a pore  of  dia- 
meter d is  'given  by 


(1-g) 

for  the  1st  layer 


f2 

N 


Tl/g 


10 

2TT075 


./N 


(4) 

NMR  measurements  have  shown  that 
-7 


T 


(//  2dVp  + A2  T / 2d2  ) , 


(20) 
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t he  first  term  allowing  for  the  cosine-law  of  desorption  then  second  for  the  delay 
which  is  related  to  the  sticking  time  T . The  most  general  expression  for  diffusion 
in  the  regime  A >>  d,  including  the  effects  of  surface  diffusion,  is  then 


= “ fr*2  + 3td8) 

3 (d/VF  + j ) 


(21) 


where  Ds  is  a surface  diffusion  term.  In  this  first  approximation  we  Ignore  the  effects 
of  surface  diffusion  in  the  second  layer  (setting  D and  r = 0)  and  take  the  second  layer  con- 
tribution  D0  to  be  given  by  equation  (19)  with  V = o,  namely 

°2  " i <’  + V)dVF  <22> 


the  Z term  allowing  for  the  fact  that  we  have  a quasiparticle  system.  The  first  layer 
contribution  is  given  by 


(23) 

2 . 2 


Di  ■ i d2L1  + f f h y<  £ + t> 

c d F 

“7  1 1 r” 

We  expect  t ^ T„  ~ 1 0 secs,  so  that  r » d /V  , giving  D = * — + 7 — 

^ C f 1 a x ^ T 

1 C 

« D, Therefore  spin  diffusion  is  not  important  in  the  first  layer 


w hich  is 


(D1  = 0). 


f)  The  net  diffusion  coefficent  is  then  given  by 

D " aiD1  + a2°2  + af°f 
where  a1  = ci  xi  + °(B) 

a2  = fX2(C2  + Cf)/X  8-  0 (24) 

af  = (c2  X 2 + Cf  Xf>  ' * 

with  ai  + a2+af  = 1 (defining  X as  In  the  two-fluid  susceptibility  model 

(section  2).  In  the  expressions  for  a^  we  have  set  g = 0 for  simplicity  so  that  the 
term  is  <<  the  other  terms  and  the  sole  function  of  the  first  layer  is  in  the  weighting 
of  the  overall  expression  by  \ Obviously  g cannot  be  identically  zero  if  the  model  1b  to 
be  self  consistent.  One  is  considering  (in  principle)  exchanges  with  the  first  layer. 

The  definition  of  &2  is  straight-forward.  The  definition  of  a^  allows  for  particles  leav- 
ing the  second  layer  and  diffusing  in  the  central  region,  by  bulk  processes,  and  then 
returning  to  that  layer.  A further  modification  without  this  last  term  has  also  been 
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carried  out  but  makes  little  difference  to  the  final  result.  In  both  models  we  have  assumed 

3 

Dj.  equals  the  diffusion  coefficient  of  bulk  liquid  lie  . The  curves  show  that  good 
agreement  is  obtained  with  a scaling  factor  a 0. 15,  consistent  with  e = . 3^  x ~ .5 
(figs^).  We  note  that  a slight  maximum  is  predicted  around  0.  07  K which  is  reproduced 
in  the  results.  However,  the  value  of  D is  also  predicted  to  decrease  at  low  T which  does 
not  appear  to  be  consistent  with  the  data.  The  same  analysis  could  be  carried  out  with 
non  zero  values  of  g.  although  on  physical  grounds  one  expects  g to  be  small;  this 
could  however  be  important  in  that  the  main  reason  for  the  drop  in  D predicted  is  the 
rise  in  c x • in  relation  to  the  total  x*  without  a compensating  term  a^  D^. 

Qualitative  predictions,  which  might  be  used  to  test  the  above  model  experi- 
mentally are:- 

(i)  a drop  in  spin  diffusion  coefficient  at  lower  temperatures. 

4 

(ii)  a considerable  increase  in  D if  the  first  layer  is  replaced  by  He,  but  with 
no  significant  drop  appearing  at  lower  temperatures.  An  increase  might 
also  be  expected,  of  course,  if  the  specular  reflection  coefficient  increase^ 
but  comparisons  between  D and  mass  flow  experiments  should  allow  dis- 
crimination between  these  two  possibilities;  if  v changes,  mass  flow  and  D 

should  be  affected  to  the  same  extort#;  If  the  explanation  above  is  correct,  mass 

(52) 

flow  will  not  be  significantly  changed'  . 

Spin  Diffusion  Allowing  for  Surface  Diffusion 


The  previous  analysis  of  the  spin  diffusion  experiments  ignored  the  possibility 
of  quasiparticles  diffusing  in  a two-dimensional  sense  within  ,the  second  layer.  While 
it  is  true  of  course  that  the  quasiparticle  concept  is  probably  not  valid  in  such  situations 
it  seems  worthwhile  for  two  reasons  to  examine  the  possibility  phenomenologically. 

First,  in  gas  experiments  at  high  temperatures,  it  is  usually  the  case  that  surface  dif- 
fusion dominates  Knudsen  diffusion.  Second,  if  we  assume  that  the  surface  diffusion 

term,  D , is  similar  to  the  bulk  diffusion  coefficient  at  30  atm,  then  a large  con- 
8 

tribution  will  be  expected  at  low  temperatures.  The  effect  of  D will  over-corn pensate 
for  the  drop  due  to  the  changing  susceptibilities  of  the  layers,  causing  D to  eventually  rise. 
Such  an  effect  could  be  looked  for  experimentally. 


The  expression  for  D_  is 

V 

fl/3  (1  ♦ Z„/4)  Vp  d * CJ-)  D,  ] 
°2  ’ [1  ♦ VF^d] 


where  1 + Z /4  = . 33  at  S.  V.  P.  ( j - 5 x 10  sec;  V 

o * 


= 5. 38  x 10  ctn/sec; 


- 39  - 


d = 60$)  and  D = D (30  atm)^\  The  total  D is  then  given  by  equation  24. 
The  resulting  curve  (fig  ^ ) goes  through  a maximum  at  07  K ^or  reasons  men- 
tioned earlier)  but  instead  of  decreasing  monatonically  continues  through  a mini- 
mum at  ~ 10  - 15  mK  and  increases  thereafter. 

At  very  low  temperatures  equation  (25)  simplifies  to 

D - VC2  + CU  X2/X 

where  x 2?  c-j  and  XQ  = 2.95T^8\ 
giving 

D = aT  + b/T  (26) 

with 

a = 1/3  (1  + Zq/4)  VFd(c2+cL)(X2/X)/(l  + VFr/d) 

b = (VFr/d)D8  (c2  + cL)(X2/X)/(l  + VFr/d) 


A minimum  is  expected  when  T 


= (b/a)i 


T 


min 


f 3 TPs 

L 2 , 

d (1+Zq/4) 


4/2 

J 


(27) 


Using  c = 0.28,  c = 0.18  and  c = 0.54^|t  = 3x10  **  seconds  from  equation  3 

‘ 2 L*  g 

and  assuming  that  Dg  is  equal  to  the  spin  diffusion  coefficient  of  bulk  He  at  ~30  atm.  It 

Ls  found  that  T . ~ 10  mK.  A smaller  mean  pore-diameter,  or  considerable  number 
min  ~ 

of  smaller  pores,  will  raise  T^^.  The  behaviour  predicted  by  equation  (26)  would  of 
course  be  open  to  experimental  study. 
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flow  measurements  (if  the  gas  accommodation  coefficient  is  i).  (ill)  in 
the  mass  flow  experiment  particles  travel  over  macroscopic  distances  and 
therefore  sample  the  whole  spectrum  of  bends,  constrictions  and  blind  holes. 
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29.  In  the  thermal  conductivity  experiment,  there  is  the  possibility  of  thermal 
short  circuits  through  the  Vycor(31),  arid  in  the  X « d regime  thermal 
energy  is  passed  on  by  collisions  so  that  the  particles  themselves  need  to 
travel  relatively  small  distances.  If  in  the  spin  diffusion  experiment  the 
characteristic  distance  (which  is  related  to  the  size  of  the  field-gradient'  ) 
is  comparable  with  d,  then  X will  be  closer  to  unity.  Since  D increases  some- 
what at  lower  temperatures,  the  scaling  factor  could  even  be  regime-dependent. 


30. 


Expressing  T 1 in  Kelvins,  we  find  that  the  rate  of  Interaction  is  ~1K, 
comparable  with  the  value  for  bulk  3He  at  S.  V.  P. 


31.  I am  grateful  for  conversations  with  A,  J.  Leggett,  M.  Moore  and  Shin 

Takagi,  all  of  whom  however,  have  expressed  some  considerable  dissent'. 
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Phys.  Letters  47A,  215  (1974). 
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Physics  Letters. 

38.  See,  for  example,  the  work  of  K.  Schiffman. 

39.  A.J.  Leggett,  private  communication. 

40.  A.  Eggington  and  M.  M.  Moore, 
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The  authors  have  used  a similar  idea  to  explain  the  enhanced  specific 
heat  of  a dilute  3IIe  in  4He  mixture  contained  within  Vycor  pores. 

41.  M.G.  Richards  has  used  similar  ideas  for  explaining  Tj  data  in  gaseous 
3He:  M.G.  Richards  and  R,  Chapman,  Phys.  Rev.  Lett.  33,  18  (1974). 

42.  The  melting  pressure  drops  very  slightly  as  field  increases  (W.  Truscott, 
private  communication). 

3 4 

43.  The  initial  recovery  for  a second  layer  of  He  absorbed  on  He  is  believed 
to  occur  with  a time  constant  250  msecs  at  6 MHa(^).  Since  this  layer  is 
expected  to  be  very  mobile,  Tj  is  likely  to  be  Independent  of  frequency.  If 
the  second  layer  remains  Intact  when  more  layers  are  present,  then  T.  for 
the  saturated  film  may  level  off  for  frequencies  greater  than  6 MHZ.  on  the 
other  hand,  if  that  layer  becomes  more  liquid-like  (long  Tj),  then  the  mono- 
layer  will  contlne  to  be  the  relaxation  zone. 
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44.  H.A.  Reich, 

Phys.  Rev.  129,  630  (1963). 
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45.  Similar  effects  occur  in  solid  He-  He  mixtures;  a review  by  M,G,  Richards 
is  in  publication. 

46.  The  solution  of  the  diffusion  equation,  with  relaxation  occurring  at  the  walls, 
also  leads  to  non-exponential  recovery^), 

47.  This  situation  is  similar  to  the  gaseous  counter  diffusion  experiments  of 
Basmadjian  and  Chu  (26), 

48.  (i)  Electron  micrographs  of  VycoTshow  a factor  of  two  variations  in  pore 
diameter,  K.  Atkins  et  al  Phys.  Rev.  102,  582  (1966). 

(ii)  In  reference  26,  the  authors  propose  a model  of  cylindrical  pores, 
radius  R,  with  spherical  enlargements,  radius  2R. 

o 

(iii)  Analysis  of  desorption  isotherms  indicates  that  pores  of  radius  19  A 
must  be  present  (the  average  values  being  30A),  D.  J.  Creswell,  thesis 
^(unpublished)  University  of  Sussex. 

49.  A similar  analysis,  applied  to  the  thermal  conductivity  experiment  in  Vycor, 
gives  a very  good  fit  to  the  data  up  to  60  mk,  where  the  'shunt  conductance' 
becomes  sizeable.  The  same  value  of  F was  used  (F  = 0.09)  and  d^.  was  taken 
to  be  half  the  volumetric  diameter  (dc  = 30A). 

50.  With  a slightly  reduced  scaling  factor  it  should  be  possible  to  use  the  same 
model  with  dp  = 70A  (Instead  of  dg  = 36A),  If  the  scaling  factor  is  chosen  to 
give  agreement  at  low  temperatures,  the  data  for  0. 1K-0.2K  would  lie  above 
the  bulk  liquid  curve  by  an  amount  less  than  the  experimental  error  («60%). 
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if  the  first  3He  layer  is  replaced  wi  th  4He. 

53.  A.C.  Anderson,  W.  Reese  and  J.C.  Wheatley, 

Phys.  Rev.  127,  671  (1962). 

54.  Experiments  carried  out  by  J.  Holt,  M.G.  Richards,  D.  J.  Creswell  and 
A.L.  Thomson  (June  1972). 

55.  I am  grateful  to  Alias  Jane  Thornton  for  her  help  in  computation  and  graph- 
plotting. 
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Fig.  2.  Initial  spin-lattice  relocation  lime  versus 
frequency. 
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Fig.  3.  T ^ plotted  against  (frequency)  for  a saturated  He  film  (full- 
pores).  The  temperature  was  approximately  0. 1 K.  (reference  54). 

T Milliseconds 

N w w 00  vo 

Oooooo  oo  oo 


FREQUENCY  (MHz) 


*> 
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approximately  0.1  K.  (reference  54). 
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Fig.  6.  T as  n function  of  temperature  for  a full-pores  film  (from  reference  10). 
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*'lg.  ft  The  spin  diffusion  coefficient  of  He  in  Vycor,  compared  with  a model 
which  allows  for  delaying  effects^®). 
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Fig.  ^ I'ho  spin  diffusion  coefficient  of  ’ He  in  Vycor,  compared  with  a model 
allowing  for  surface  diffusion^- . 
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□ I.  Minima  in  Spin-Lattice  Times  in  Adsorbed  Helium 
A.  L.  Thomson,  D.  F.  Brewer  and  Y.  Goto. 

University  of  Sussex,  School  of  Mathematical  and  Physical  Sciences, 
Kalmer,  Brighton  BN1  9QI1,  England. 


In  previous  papers1  we  have  described  several  measurements  of  the  spin- 

lattice  (T  ) and  spin-spin  (T  ) relaxation  times  of  *IIe  films  adsorbed  on  Vvcor 

porous  glass  within  the  temperature  range  0.3  T < 1.5K,  at  various  coverages 

up  to  one  monolayer.  The  ai  m of  the  present  experiments  was  to  explore  these 

more  fully  over  a more  extensive  range  of  frequency  (0.50  M 1 1^  to  2.00  M1M 

an/1  temperature  (0.3K  to  4.2Kb  We  have  found  minima  in  the  spin-lattice 

recovery  time  as  a function  of  temperature  which  are  similar  to  those  occurring 
2 3 

in  bulk  solid  ’ . 

Details  of  the  apparatus  have  been  reported  elsewhere1  and  therefore  only 

3 

the  briefest  outline  is  given  here.  The  sample  helium  amounting  to  32  cm  NTP 

. 2 , 

and  thought  from  adsorption  isotherm  measurements  to  correspond  to  complete 
monolayer  was  adsorbed  on  Vycor  glass  contained  in  a container  in  a chamber  made 
of  Stycast.  Around  this  chamber  was  wound  a copper  coil  whose  resonant  frequency 
was  590KH  and  higher  frequencies  were  obtained  by  connecting  a parallel  with  this 

7. 

coil  a high  Q coil  located  outside  the  cryostat.  A pulsed  nuclear  magnetic 

resonance  technique  was  used  to  measure  the  relaxation  times,  employing  a 

simple  180°-90"  pulse  sequence  for  and  a 90°~180'’  sequence  for  T^.  Echoes 

were  obtained  corresponding  to  different  separation  times  between  the  pulses  and 

the  heights  of  these  echoes  were  recorded.  The  echoes  observed  with  the 

measurements  always  showed  a good  exponential  recovery  but  this  was  not  the 

case  for  the  measurements  where  the  recovery  had  a characteristic  time 

which  lengthened  as  the  signal  moved  towards  equilibrium.  Nevertheless  most 

of  the  recovery  occurred  before  this  departure  from  strictly  exponential  behaviour 

and  It  is  the  time  corresponding  to  this  recovery  that  we  take  here  to  be  Tj.  Our 

results  are  shown  in  Fig.l  where  T and  T are  plotted  against  inverse  temperature. 

1 

The  Tj  data  show  clear  minima  for  each  of  these  frequencies,  at  a temperature 
which  increases  very  slightly  with  frequency.  Below  IK,  T^  and  T^  are  both 
approximately  constant.  Above  this  temperature  Tg  increases  monotonically  up 
to  4.2K,  and  throughout  the  whole  temperature  range  T2  Is  approximately  an  order 
of  magnitude  smaller  than  Tj . A feature  of  the  Tj  data  is  that  they  increase  with 
frequency  over  the  whole  temperature  range,  while  the  T2  data  3how  no  measurable 
dependence  on  frequency. 


2 3 

There  is  a broad  similarity  between  the  data  presented  here  and  those  ’ 

for  bulk  solid  *Me.  In  the  latter  case,  the  high  temperature  region  corresponds 

to  direct  relaxation  of  the  Zeeman  energy  to  the  lattice  and  the  minimum  in  T , 

which  is  a standard  feature  of  BPP  theory,  occurs  where  the  diffusion  time  is 

around  o'  1 where  o>  is  the  Larmor  frequency.  In  the  low  temperature  plateau 
o o 

region,  the  recovery  is  dominated  by  the  Zeeman-exchange  relaxation  time,  due 
to  tunnel  ling  or  exchange  motion.  Such  a simple  interpretation  is  unlikely  to  be 
directly  applicable  to  our  observations,  since  they  are  complicated  by  the 
heterogeneous  nature  of  the  glass  substrate  which  may  account  for  the  non-exponential 
character  of  the  spin-lattice  recoveries.  However,  it  is  worthwhile  examining  the 
implications  for  diffusion  limes  and  activation  energies  in  the  high  temperature 
region,  and  for  the  characteristic  tunnelling  frequency  in  the  low  temperature 
region. 

VVc  note  that  the  temperatures  of  the  minima  increase  with  increasing 
frequency  and  occur  at  1 . 92  i.  02K,  1.96  ± . 02K  and  2.00  ±.02K.  If  we  ascribe 
a diffusive  time  to  the  motion  and  an  energy  of  activation,  such  that 
T - Tq  exp  («  /r  ) 

then  we  can  by  fitting  these  values  obtain  values  of  tq  and  The  activation  energy  woulc 

appear  to  be  of  order  30K  while  the  constant  r ha9  a value  of  order  10  ***  secs. 

o 

Such  a large  value  for  the  activation  energy  does  not  seem  to  be  realistic,  however, 

since  it  would  imply  a much  sharper  minimum  than  observed.  To  i nterpret  the 

plateau  region,  one  requires  knowledge  of  the  form  ot  the  spectral  density  describing 

the  motion  of  the  spins.  Karlicr  measurements1  suggested  that  a Lorentzian 

spectral  density  (corresponding  to  an  exponential  correlation  function)  might  be  the 

most  appropriate,  and  certainly  the  comparatively  small  variation  of  T^  with 

frequency  is  in  accord  with  this.  An  analysis  of  these  data  is  currently  being  made 

assuming  various  spectral  densities  but  an  estimate  for  the  time  characterising 

the  proposed  tunnelling  motion  can  be  found  from  the  relationship 

T - l/T-Ao* 
tun  2 

which  holds  for  a Lorentzian  spectral  density  in  the  region  where 

, r~~2 

uT 
tun 

is  the  second  moment  of  the  resonance  line,  and  can  be  calculated  if  the  relative 
positions  of  neighbouring  spins  are  known.  One  extreme  for  A o'  would  occur  if 
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the  spins  were  uniformly  distributed  which  for  our  adsorbate  would  produce  a 
value  of  'lO^sec  “.  Another  extreme  would  be  tight  packing  of  the  atoms  into 
islands  composed  of  a triangular  lattice  such  that  the  atoms  were  at  the  spacing  of 
the  minimum  of  their  van  der  Waals  potential.  This  would  produce  a value  for 
£ a'2  of  ~ 5x1 0^  sec  2.  Thus  r would  be  expected  to  lie  within  the  range  10 
to  10">  secs.  Further  measurements  should  elucidate  this  picture  and  in 
particular  experiments  are  planned  for  several  diffei'ent  coverages  of  adsorbate. 

Supported  by  SRC  (Grant  B/RG/1795)  and  US  Army  (DA-ERO-1 24- 74- (#0046). 

References: 

1.  D.  F.  Brewer,  D.  J.  Creswell,  Y.  Goto,  M.  G.  Richards,  J.  Rolt  and 
A.  L.  Thomson.  Monolayer  and  Submonolayer  Helium  Films,  Ed. 

J.  G.  Daunt  and  E.  Lerner  (Plenum  Press  1973)  p.101  and  references 
therein. 

2.  R.  C.  Richardson,  E.  Hunt,  H.  Meyer.  Phys.  Rev.  m A1326  (1965) 

3.  M.  G.  Richards,  J.  Hatton  and  R.  P.  Gifford,  Phys.  Rev. 139  A91  (1965) 


